Single-qubit lasing and cooling at the Rabi frequency 
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For a superconducting qubit driven to perform Rabi oscillations and coupled to a slow electro- 
magnetic or nano-meclianical oscillator we describe previously unexplored quantum optics effects. 
When the Rabi frequency is tuned to resonance with the oscillator the latter can be driven far from 
equilibrium. Blue detuned driving leads to a population inversion in the qubit and a bi-stability 
with lasing behavior of the oscillator; for red detuning the qubit cools the oscillator. This behavior 
persists at the symmetry point where the qubit-oscillator coupling is quadratic and decoherence 
effects are minimized. There the system realizes a "single-atom-two-photon laser". 
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Several recent experiments on quantum state engineer- 
ing with superconducting circuits P, [3, 0] realized con- 
cepts originally introduced in the field of quantum optics 
and stimulated substantial theoretical activities 0, 
Josephson qubits play the role of two-level atoms, while 
oscillators of various kinds replace the quantized light 
field. Motivated by one such experiment [l|, we inves- 
tigate a Josephson qubit coupled to a slow LC oscilla- 
tor (Fig. [T] a) with eigenfrequency (in the MHz range) 
much lower than the qubit 's energy splitting (in the GHz 
range), lot ^ A.E. The qubit is ac-driven to perform 
Rabi oscillations, and the Rabi frequency is tuned 
close to resonance with the oscillator. For this previously 
unexplored regime of frequencies we study both one- 
photon (for Qjf w lot) and two-photon (for flji w 2lot) 
qubit-oscillator couplings. The latter is dominant at the 
"sweet" point of the qubit, where due to symmetry the 
linear coupling to the noise sources is tuned to zero and 
dephasing effects are minimized ^ . When the qubit driv- 
ing frequency is blue detuned, Slo = LUd ~ AE > 0, we 
find that the system exhibits lasing behavior; for red de- 
tuning the qubit cools the oscillator. Similar behavior 
is expected in an accessible range of parameters for a 
Josephson qubit coupled to a nano-mechanical oscillator 
(Fig. [T]d), thus providing a reahzation of a SASER 7] 
(Sound Amplifier by Stimulated Emission of Radiation). 

Our work is motivated by a recent experiment [l| where 
under the (assumed) condition flu ^ lot a strong excita- 
tion of the LC oscillator was observed. Previous attempts 
to explain this behavior did not resolve several is- 

sues: First, according to these theories the effect should 
vanish at the symmetry point, but it was observed there 
as well. A possible explanation P, Q , assuming uncon- 
trolled small deviations from the symmetry point, is not 
supported by experiments. Here we explore an alterna- 
tive, namely that 2-photon processes at flu « 2u;t which 
persist at the symmetry point, are responsible for the 
observed effect. 

The second unresolved problem is the magnitude of the 
effect. The experiment showed an increase by a factor 
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FIG. 1: The systems, a) In the circuit QED setup of Ref. [ll] 
an externally driven three-junction flux qubit is coupled in- 
ductively to an LC oscillator, b) In an equivalent setup a 
charge qubit is coupled to a mechanical resonator. 



~ 10 in the number of oscillator quanta [H. The theory 
of Rcf. [8], valid in the perturbative regime, predicts a 
much weaker effect. Here we demonstrate that the high- 
frequency driving of the qubit, in combination with the 
unusual resonance condition introduces qualitatively new 
effects: Namely, for blue detuning of the qubit driving a 
population inversion is created between the qubit dressed 
which are split by the Rabi frequency [lol . 



states, which are split by the Rabi frequency [10|, lll|. 
Under these conditions the ac-driving provides more en- 
ergy than needed to excite the qubit, and excess energy 
is fiowing into the resonator. The lasing threshold and a 
proper lasing state with the characteristic line narrowing 
can be reached, and the system becomes a "single-atom- 
laser" fl3] . Similarly, for red detuning we predict strong 
cooling of the oscillator. 

The described mechanism is strong enough to explain 
the observed level of excitation of the oscillator. At this 
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stage, using parameters provided by Ref . [l| , we have not 
reached quantitative agreement with experiment, possi- 
bly because of a combination of strong drive and low 
quality factor. However, the effects should be observable 
in circuit QED setups of the type described if sufficiently 
high Q-factor can be achieved. Furthermore, our analy- 
sis applies for low-frequency nano-mechanical resonators, 
where sufficiently high Q- factor have already been ob- 
tained, and true lasing (sasing) should arise. 

The systems considered are shown in Fig. [TJ To be 
specific we first analyze the Rabi driven flux qubit cou- 
pled to an LC-oscillator (Fig. [T^) with Hamiltonian 

= - i e {^f) a^-^Aaa:- f^i?o cos(wdi) 

+ + 9<7z{o' + a^) ■ (1) 

The first two terms describe the qubit, with Pauli ma- 
trices acting on the flux basis states of the qubit. The 
energy bias e(<&^'^) between the flux states is controlled 
by a dc magnetic flux, and A is the tunneling amplitude 
between both. The third term accounts for the driving of 
the qubit by an applied ac flux with amplitude ^Irq and 
frequency tOd- The last two terms describe the oscillator 
with frequency ujt = {LtCt)~'^^'^ , as well as the qubit- 
oscillator interaction. We estimate the coupling constant 
as g ~ MIplT.Oi where M is the mutual inductance, Ip 
the magnitude of the persistent current in the qubit, and 
It,o = [y^T I'^Tt)^!'^ the amplitude of the vacuum fluctu- 
ation of the current in the LC-oscillator. 

After transformation to the eigenbasis of the qubit, 
a Schrieffer- Wolff transformation, and a rotating wave 
approximation (RWA) the Hamiltonian reads 

H = AE (Tz + cos (uJdt) cos (Tx + ujt o,^ ci 

9 
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+ gsinCcFz {a + a^) - -^cos^ Cc^z {a + a'') (2) 

where AE = ve^ + A^ and tanC = e/A. The derivation 
resembles that of Ref. [H] in the dispersive regime, except 
that we keep terms proportional to and . They 
cannot be dropped, since the oscillator's frequency, ujr, 
is much lower than AE. Hence, and evolve slowly. 

Next we transform the Hamiltonian ([2]) into the rotat- 
ing frame, which introduces the full Rabi frequency ^Ir — 
\/^j^Q cos^ C + ^"^^i and an angle tan/3 = (5a;/(riflo cos^) 
characterizing the detuning Suj = tOd — AE. We employ 
a second RWA in which we consider lot and Qr ~ lot as 
fast (see justification below). In the interaction represen- 
tation with respect to the non-interacting Hamiltonian 
Ha = {VIr/2) (Jz + LOT a^a^ we arrive at 

-I- 53 {a) a + aa^) Oz . (3) 



We kept both single- and two-photon interactions with 
5i = gsin(^cos/3 and (72 = (5^/A) cos^ C cos/3, although 
within the RWA only one of them survives: the first one 
near the resonance at fi/f ~ lot^ the second one for fJ/j 
2lot- The third term with 53 = — (q^ / A) cos^ C sin (3 
causes a frequency shift of the oscillator [13| . All three 
coupling constants are much smaller than the bare g in 
the vicinity of the symmetry point. In what follows we 
assume that the qubit is kept near the symmetry point, 
i.e., e ^ A and cos^ ~ 1. Exactly at the symmetry point 
the single-photon coupling constant gi vanishes, but the 
two-photon couplings 52 and 53 persist. 

To account for the effects of dissipation we in- 
troduce two damping terms in the Liouville equa- 
tion for the density operator of the system, p = 
— i [H,p\ -\- Lq p + Lup. As far as the qubit is con- 
cerned we consider only spontaneous emission with rate 
Fq. Hence Lq p = Fq {2a^pa^ — pa+o- — cr+a-p) /2. 
This restriction is justified as long as the tempera- 
ture is lower than the qubit's energy splitting, and 
the system is biased near the degeneracy point (e = 
0) where the additional "pure" dephasing is weak. 
The resonator damping can be expressed in standard 
form Q, ,0 = K (A^th + 1) {2apa^ - a^ap - pa^ a) /2 + 
KNth {2a^ pa — aa^ p — paa^) /2, where k denotes the res- 
onator width and Nth = 1/ [exp(wT/fc_BT') — 1] the ther- 
mal average number of photons in the resonator. Since 
we consider qubits with high degree of coherence, we are 
motivated to study the regime Fq, k, 51, 32, <?3 ^ <^t, ^r, 
where the second RWA performed above with respect to 
frequencies lot, ^r is justified. 

After transformations to the rotating frame and the 
interaction representation, within the RWA the qubit 
damping takes the form 

Lqpf = ^{2a-pf a+- pf<J+a--a+a-pf ) 

+ ^-^[azP^az-p"}) . (4) 

Remarkably, although we started from spontaneous emis- 
sion only, the transformation to the rotating frame in- 
troduces transition in either direction with rates F|_| = 
^(1 ± sin/3)^ and an additional pure dephasing term, 
F* = cos^ /3. Right on resonance, where /3 = 0, we 
have F| = Fj^, corresponding to infinite temperature. We 
further note that detuning modifies the effective temper- 
ature in the rotating frame. For blue detuning, /? > 0, 
we find F| > F^, i.e., negative temperature. 

The Hamiltonian ([3]) combined with the Liouville 
equation and qubit damping term (|4]) provides the 
the theoretical description of the system. It differs 
from earlier work in the unusual resonance condition 
with the Rabi frequency and the inclusion of both 
1- and 2-photon processes. The latter dominate at 
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the symmetry point, where decoherence effects are 
minimized. Below we analyze both resonance situ- 
ations Ofl LOT or rjfl '^ojt, and investigate the 
effects of blue or red detuning, 5uj = lo^ — AE, of the 
qubit driving frequency. We proceed in the frame of the 
standard semi-classical approach IJ, ll5| of laser physics. 



(i) One-photon interaction. When the Rabi frequency is 
near the resonance with the oscillator, fi^; « lot, we can 
neglect in RWA in the Hamiltonian ([3]) the term with 
coupling constant g2- If we - to start with - neglect 
fluctuations, the system is described by Maxwell-Bloch 
equations for the classical variables a = {a), a* = {a^), 
s± = ((T±), and Sz = (o'z), which can be derived from 
the Hamiltonian if all correlation functions are assumed 
to factorize. The qubit variables can be adiabatically 
eliminated as long as k, gi ^ Fq, which leads to a closed 
equation of motion for a. If we account for fluctuations, 
e.g., due to thermal noise in the resonator, a becomes a 
stochastic variable obeying a Langevin equation [3| , 
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2 + (5) 



Here C = 2gj, sf = -£)o/ (l + |a^|/"o) is the sta- 
tionary value of the population difference between the 
qubit levels, and Dq = {Ti — T-\)/Ti is the normal- 
ized difference between the rates with Fi = F| + F|. 
We further introduced the photon saturation number 



no = T^pVi/Agi and fio 



= no(l + (5f2Vr^), and the 
total dephasing rate F^^ = Fi/2 + F* . The detuning 
of the Rabi frequency enters in combination with a fre- 
quency renormalization, 50, = Q.^ — i^r + 4^3 [a p. The 
Langevin force due to thermal noise in the oscillator satis- 
fies {C{t)C{t')) = KNth6it-t') and = 0- Noise 
originating from the qubit can be neglected provided the 
thermal noise is strong, KiVth ^ gl/T^. 
(a) Two-photon interaction. The two-photon effect dom- 
inates near the resonance condition « 2lot- In this 
case in RWA we neglect in ^ the term with coupling 
constant gi. The corresponding Langevin equation for 
the resonator variable is of the same form as Eq. ^ 



but with C = 4g||a| and 



-Do/(l + {\a^\/rior 



The photon saturation number is now given by uq ~ 
(F^Fi/4(7|)V2, and no = no(l + SnyT^^/^ Again 
^(i) represents thermal noise, while noise arising from 
the qubit can be neglected if KiVth ^ gln/F,^, where 
n = (|ap) is the average number of photons. The rele- 
vant detuning of the Rabi frequency for two-photon in- 
teraction is given by SV, = flff, — ^lut + 4^3 |ap. 

In Fig. [2] we summarize our main results obtained by 
solving the Langevin (Fokker-Plank) equations 15| . The 
number of photons n is plotted as a function of the de- 
tuning Slu of the driving frequency and driving ampli- 
tude r^ijo- It exhibits sharp extrema along two curves 
corresponding to the one- and two-photon resonances. 
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FIG. 2: Average number of photons in the resonator as func- 
tion of the driving detuning 5lj and amplitude S7flo. Peaks 
at 5cij > correspond to lasing, dips at 5lj < to cool- 
ing. The inner curve is the one-photon resonance, which 
exists only away from the symmetry point. (Here we as- 
sume e = O.OIA.) The outer curve is the two-photon res- 
onance, which persists at e = 0. In domains of bi-stabihty 
the lowest value of n is plotted. Sharp drops in the curves 
are related to the bi-stability bifurcations. The parame- 
ters for the qubit are A/27r = 1 GHz, e = O.OIA, and 
ro/27r = 125 kHz, the frequency and line-width of the res- 
onator are ut/2-k = 6MHz and Hi/2n = 1.7 kHz, the coupling 
constant is g/2Ti — 3.3 MHz and the temperature of the res- 
onator r = 10 mK. The inset shows the bistability of the 
photon number for Q,ro/2tt = 7MHz. The dashed line repre- 
sents the unstable solution. 



^B. = — 4(7371 and — 2ujt — 4g3n. Blue detun- 
ing, 5lj > 0, induces a strong population inversion of the 
qubit levels, which in resonance leads to one-qubit lasing. 
In experiments the effect can be measured as a strong 
increase of the photon number in the resonator above 
the thermal values. On the other hand, red detuning 
produces a one-qubit cooler with photon numbers sub- 
stantially below the thermal value. Near the resonances 
we find regions of bi-stability illustrated in the inset of 
Fig. [2] In these regions we expect a telegraph-like noise 
due to random switching between the two solutions. 

For strong population inversion, F| <C F|, i.e., Dq ~ 
— 1, we obtain h ~ iVth -I- Fi/(2k) for single-photon 
and n ^ A^th + Fi/k for two-photon resonances. If 
Ti/k ^ A'th we obtain the usual lasing state with Pois- 
son distribution p„ and corresponding line narrowing. In 
the opposite limit the state is almost thermal. 

For comparison we plot our results using parameters 
as provided in Ref. As shown in Fig. [5] blue detuning 
causes a increase in the resonator occupation number by 
a factor ~ 3, a value which is short of the experimentally 
observed one. By adjusting the parameters we can reach 
enhancements comparable to the observed ones. It re- 
quires either decreasing k by a factor of 4 or increasing 
Fq by a factor of ~ 4 or a combination thereof. Further 
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analysis and fitting as well as studies of higher-Q circuits 
should help resolving the remaining discrepancies. 

So far we described a flux qubit coupled to an LC 
oscillator, but our analysis applies equally to a nano- 
mechanical resonator capacitively coupled to a Joseph- 
son charge qubit (see Fig. [1)3). In this case az stands 
for the charge of the qubit, and both the couphng to 
the oscillator and the driving are capacitive, i.e., involve 
Uz- To produce capacitive coupling between the qubit 
and the oscillator, the latter is metal coated and charged 
by a voltage source The dc component of the gate 
voltage Vg puts the system near the charge degeneracy 
point where the dephasing due to the 1// charge noise 
is minimal. Rabi driving is induced by an ac component 
of Vg. Realistic experimental parameters are expected 
to be very similar to the ones used in the examples dis- 
cussed above, except that a much higher quality factor 
of the resonator (~ 10^) and a much higher number of 
quanta in the oscillator can be reached. This number will 
easily exceed the thermal one, thus a proper lasing state 
with Poisson statistics, appropriately named SASER 
is produced. One should then observe the usual line nar- 
rowing with line width given by KA'th/(4n) ~ K^TYth/Fi. 
Experimental observation of this line-width narrowing 
would constitute a confirmation of the lasing/sasing. 

Another useful application of the considered scheme is 
the cooling of the resonator. From our analysis we con- 
clude that a population of order n — 1 can be reached 
for optimal detuning. Further analysis is required to de- 
termine the optimum conditions. 

The effect which we described here differs from the 
"dressed-state lasing" studied earlier in quantum op- 
tics where the resonator is coupled to the high- 
frequency MoUow fill transitions with energy differences 
/S.E ± ilfl. In the present setup the oscillator couples 
to the much lower Rabi frequency (similar to Ref. [iq). 
which can be readily tuned to resonance with the oscil- 
lator in order to reach the lasing threshold and a proper 
lasing state. 

Laser-like behavior has been predicted also for nano- 
mechanical resonators coupled to superconducting single- 
electron transistors [l7] . Our proposal differs in two ma- 
jor points: (i) As an active medium we have a coherently 
driven qubit instead of a dissipative transistor. While 
transistors usually provide a broad band negative tem- 
perature environment for the oscillator, frequencies of or- 
der LOT are resolved in our situation. Experience gained 
in the field of quantum optics, e.g. from the comparison 
of Doppler cooling and resolved side-band cooling, shows 
that frequency resolved regimes lead to stronger effects. 



(ii) We consider two-photon coupling which allows las- 
ing at the degeneracy point of the qubit where the sen- 
sitivity to low-frequency noise is substantially reduced. 
These facts provide sufficient flexibility in the choice of 
system parameters and should allow reaching a proper 
lasing state with photon numbers much above the ther- 
mal value and the characteristic line- width narrowing. 
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